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ABSTRACT 

The role of thermal conduction in regulating the thermal behavior of cooling flows in galaxy clusters 
is reexamined. Recent investigations have shown that the anisotropic Coulomb heat flux caused by a 
magnetic field in a dilute plasma drives a dynamical instability. A long standing problem of cooling 
flow theory has been to understand how thermal conduction can offset radiative core losses without 
completely preventing them. In this Letter, we propose that magnetohydrodynamic turbulence driven 
by the heat flux instability regulates field-line insulation and drives a reverse convective thermal flux, 
both of which may mediate the stabilization of the cooling cores of hot clusters. This model suggests 
that turbulent mixing should accompany strong thermal gradients in cooling flows. This prediction 
seems to be supported by the spatial distribution of metals in the central galaxies of clusters, which 
shows a much stronger correlation with the ambient hot gas temperature gradient than with the parent 
stellar population. 

Subject headings: cooling flows — galaxies: clusters: general — magnetic fields — MHD — instabilities 
— turbulence 



1. INTRODUCTION 

Although the intracluster medium (ICM) of galaxy 
clusters has been an active area of study since the ear- 
liest days of X-ray astronomy (Felten 1966), the com- 
plexity of the physical processes regulating the behavior 
of this gas continues to surprise us. X-ray observations 
imply that cooling times in the core regions of ICM atmo- 
spheres are significantly shorter than cluster ages (Peter- 
son & Fabian 2006). Prior to the launch of XMM-Newton 
and Chandra, it was widely held that these short cool- 
ing times lead to flows in which mass is dropping out 
(Fabian 1994). Subsequent XMM-Newton spectroscopy 
confirmed the short cooling times and inwardly decreas- 
ing ICM temperature profiles, but also revealed a lack of 
gas cooling below temperatures of T ~ Tvir/S, where Tvir 
is the virial temperature of the cluster (Peterson et al. 
2001; Tamura et al. 2001). This observation, combined 
with a realization that unchecked cooling flows would 
lead to central galaxies far above the high-end cutoff of 
the luminosity function (Benson et al. 2003), made it 
clear that the cooling cores of galaxy clusters must be 
continuously heated. 

A new paradigm consisting of core heating by a cen- 
tral radio galaxy has taken hold. Although a few exam- 
ples were known at the time of ROSAT, Chandra and 
XMM-Newton have observed cavities and weak shocks in 
the ICM of many clusters that are driven by the central 
radio-galaxy (e.g., Fabian et al. 2000, 2003; Young et 
al. 2002; Heinz et al. 2002). This has motivated the hy- 
pothesis that core losses are offset by thermalization of 
the kinetic luminosity of the jets from the active galactic 
nucleus (AGN). In some (but not all; see §4) systems, 
the energetics of the AGN is comparable to the radiative 
luminosity of the cluster core (Birzan et al. 2004; Best et 
al. 2007). AGN feedback of this form has been adopted 
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by semi-analytic galaxy formation theories in the guise 
of radio-mode feedback (Croton et al. 2006). 

An alternative hypothesis is that thermal conduction 
into the core from the large heat reservoir beyond the 
cooling radius may offset the radiative loss (Binney & 
Cowie 1981). A source of uncertainty for such models is 
the factor / by which magnetic fields suppress thermal 
conduction below the classical Spitzer value ks (Spitzer 
1962). Plausible estimates range from / ~ 10~ (Chu- 
razov et al. 2001) to f ^ 0.3 (Narayan & Medvedev 
2001). Thermal conduction models have faced other 
significant challenges, however. Even unbridled Spitzer 
conductivity cannot offset radiative cooling in the lower- 
mass/cooler clusters (Voigt and Fabian 2004), a conse- 
quence of the strong temperature dependence of Spitzer 
conductivity, ks oc T^/^. For hotter clusters, thermal 
conduction can in principle offset radiative cooling if 
/ ~ 0.1 — 1 (Zakamska & Narayan 2003). In this case, 
however, fine tuning of / is required to offset radiative 
cooling, yet not erase the observed temperature gradi- 
ents within the cooling core (e.g., see Conroy & Ostriker 
2008). Furthermore, assuming / is a constant for a given 
ICM atmosphere, these equilibria are still thermally un- 
stable (Kim & Narayan 2003). 

In this Letter, we argue that recent theoretical devel- 
opments in the theory of the stability of MHD atmo- 
spheres must have a crucial bearing on the behavior of 
cluster cooling cores. In particular, it is now understood 
that in atmospheres of dilute plasmas (such as the ICM) 
in which thermal conduction is restricted to follow mag- 
netic field lines, the convective instability criterion is fun- 
damentally altered: temperature gradients, not entropy 
gradients, are crucial. Atmospheres in which the tem- 
perature gradient is decreasing upwards [dT/dz < 0), 
but field lines are isothermal, are convectively unstable 
(Balbus 2000, 2001; Parrish & Stone 2005, 2007) even 
when the entropy gradient would normally imply sta- 
bility. On the other hand, atmospheres in which the 
temperature gradient is increasing upwards {dT/dz > 0; 
the case in cooling cores), and field lines are thermally 
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conducting, are also unstable (Quataert 2008; Parrish 
& Quataert 2008). These two powerful instabilities are 
known respectively as the magnetothermal instability 
(MTI; dT/dz < 0) and the heat flux-driven buoyancy 
instability (HBI; dT/dz > 0). 

Our focus is the cooling core of a rich/hot cluster of 
galaxies. We investigate the consequences of convective 
turbulence that is driven, but also regulated, by the HBI. 
Turbulence, far from being a source of heat, is actually a 
cooling term in the energy balance. The role of convec- 
tive fluctuations is an outwardly-directed "return heat 
flux" that, together with field line insulation, lowers, but 
by no means erases, the inward diffusive thermal flux. 
This process eliminates the theoretical ad hoc fine tuning 
required of previous conduction-based models. Since our 
picture attributes the stabilization of hot cooling cores 
to the intrinsic dynamics of MHD plasmas, the observed 
homology of hot cooling cores range is perhaps not sur- 
prising. Moreover, turbulence-regulated heat conduction 
also accounts for the otherwise puzzling tight correlation 
between metalicity distribution and temperature gradi- 
ent observed in cluster cores: the two become intimately 
related to the presence of convective turbulence. 

An outline of this Letter is as follows. In §2, we discuss 
the role of radiative losses on the HBI in cluster cooling 
flows. The role of convective heat transport is studied in 
§3, and our picture for cooling cores is explicitly discussed 
in §4. We summarize our conclusions in §5. 

2. RADIATIVE COOLING EFFECTS ON THE HBI 

Before presenting our thermal model of hot cluster 
cooling cores, it is of interest to modify the HBI disper- 
sion relation of Quataert (2008) due to the presence of 
radiative losses. We restrict our calculation to the sim- 
plest and most illustrative case of a uniform background 
magnetic field in the vertical z direction. 

Radiative losses are represented by the net cooling 
function ^pC (units erg s""'^ cm~^); this terms includes 
any explicit heating processes as well as radiative cooling. 
We make the formal assumption that the background is 
static. To include first-order perturbations in radiative 
loses, the term —5{pC) should be added to the right side 
of equation (11) in Q08. Plane wave disturbances of the 
form exp(crt + ik^x + ikzz) then satisfy the dispersion 
relation 
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We have written a for —iuj in Q08, otherwise the notation 
is identical. (The quantity g is —{1/ p)dP/dz.) 

For the case of a vertical field, equation ([1]) corresponds 
precisely to the dispersion relation (13) of Q08, but with 
the additional cooling term in square brackets in equation 
([2]). The constant term in the dispersion relation ^ is 
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This should be negative for instability. In the limit of 
very small k • va, we recover the HBI condition that 
an increasing temperature profile is unstable. We see, 
however, that in the presence of a non-negligible mag- 
netic field, radiative cooling (which, in this setting, has 
{d{pC)/dT)p < 0) acts to further destabihze the atmo- 
sphere, especially at longer wavelengths. 

The non-linear behavior of an HBI unstable atmo- 
sphere is still unclear. On the basis of 2D and 3D simula- 
tions of non-cooling plane-parallel atmospheres, Parrish 
& Quataert (2008; hereafter PQ08) show that the HBI 
acts to re-orient field lines perpendicular to the tempera- 
ture gradient, effectively shutting off the conductive heat 
flux. We shall refer to this process as field-line insulation. 
However, it is not obvious that the conclusions of PQ08 
can be extended to actual cluster cooling cores. As well 
as the topological impossibility of establishing a (non- 
zero) field configuration that has zero radial field every- 
where, the inflow of matter associated with unbalanced 
radiative cooling in an insulated core will readily gener- 
ate a radial field from existing transverse field. Thus, we 
argue that field-line insulation, while surely present at 
some level, cannot be absolute in cluster cooling cores. 

One of the most salient properties of the linear HBI is 
its convective nature, in particular the fact that the prod- 
uct 6vz5T > (Q08). If sufficiently large, this reverse 
flux regulates the background heat flux: the greater the 
magnitude of dT/dz, the more vigorous the return flux. 
We suggest that the combination of field line insulation 
and reverse convective flux fixes a state of marginal sta- 
bility to the radiative HBI that corresponds to the cool 
cores of galaxy clusters. 

3. HEAT TRANSPORT WITHIN HOT CLUSTER COOLING 

CORES 

Let us now re-consider the energy equation for the cool- 
ing core of a hot cluster in the light of the above discus- 
sion. Under steady conditions, the equation of hydrody- 
namical energy conservation is 
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where $ is the potential function of the cluster, and Q = 
Ksb{b ■ VT) is the diffusive heat flux (6 is a unit vector 
in the direction of the magnetic field). 

We wish to investigate the nature of dynamically static 
solutions to this equation in which, on average, the mass 
flux pv vanishes. The HBI is assumed to be present and 
fully developed into convective turbulence. This requires 
a background heat flux Q to be present, since this is the 
seat of the instability. We assume that the density and 
temperature consist of a mean background value plus a 
fluctuation, p = pq + Sp, T = Tq + ST. The radial 
velocity v (the only component of concern here) is dom- 
inated by its fluctuating component, but must also have 
a very small second order part with a finite mean value, 
V2. This is needed if the mass flux vanishes, 

{{po + 5p){v2 + 5v)) ~ poV2 + {5p5v) = (6) 

where () indicates a time average. Defining the normal- 
ized temperature r = P/ p, our energy equation reads 



(4) 



-Po{5v 6t) 



Q 



(7) 



3 



This explicitly demonstrates that it is the divergence of 
both the (inwardly directed) diffusive heat flux Q and the 
(outwardly-directed) turbulent heat flux {5/2)p{dvST) 
that, in our model, balances radiative cooling. Notice 
that, as in Schwarzschild (1958), the convective heat flux 
is proportional to the correlation of velocity and temper- 
ature fluctuations. 

We now estimate the outwardly directed convective 
heat flux (5/2)p{SvdT), using arguments similar to 
Schwarzschild (1958). Henceforth, we drop the sub- 
script: p and T are understood to be background val- 
ues. The fastest growing modes of the HBI typically 
have (Q08) 

St ST dlnT 

— = — = $ j„ : (8) 
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where ^ is the radial displacement. Following 
Schwarzschild, we multiply this density deficit by g$,/2 
to obtain the kinetic energy per unit volume of a rising 
element, p{Sv)'^/2: 
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Therefore, the convective heat flux is 
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where P is the turbulent diffusion coefficient 
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Schwarzschild identifies ^ with 1/2, where / is a formal 
mixing length, but such details are of secondary concern 
here. The main point is that if these basic scalings are 
correct, the outward convective heat flux acts like a re- 
verse diffusive heat flux, but one in which the thermal 
conductivity is itself proportional to {dT / drY/"^ . To- 
gether with incomplete field-line insulation, this is well- 
suited to act as an efficient regulatory mechanism coun- 
teracting the otherwise large inward diffusive transport 
of thermal energy. 

4. DISCUSSION 

HBTdriven turbulence has two important conse- 
quences for heat transport in hot cluster cores, the par- 
tial suppression of the conductive heat ffux due to field- 
line re-orientation (previously stressed by PQ08) and an 
outwardly- directed convective heat flux (highlighted in 
this work). We hypothesize here that these two effects 
keep a cooling core in marginal thermal stability, with 
the (regulated) thermal conduction balancing radiative 
cooling. 

Consider a hot cluster with an initially isothermal 
ICM, lacking any embedded heat source such as an AGN. 



Cooling in the core of the cluster would set up an ini- 
tially small temperature gradient and, for a short time, 
the radiative loses would be easily compensated by ther- 
mal conduction down this temperature gradient. As time 
proceeds, however, HBI induced field-line insulation and 
the convective heat flux will act to suppress the con- 
duction heat flux, leading to a steepening temperature 
gradient. As the center of the cluster approaches a cool- 
ing catastrophe, the vigorous turbulence will comb the 
field lines of the turbulent ICM preferentially in the ra- 
dial direction, opening up an enhanced channel for the 
inwards conductive heat flux and quenching the runaway 
cooling. We posit that dense ICM cores will evolve to- 
wards a state of near static, marginal "radiative-HBI" 
stability, and identify this with the observed state of the 
hot cluster cores. 

Although this work appears to be the first to invoke 
turbulence as a cooling mechanism, the suggestion that 
cluster cooling cores may be turbulent is certainly not 
new (e.g., see Voigt & Fabian 2004). Indeed, the pres- 
ence of turbulence is motivated by observational studies 
of the rotation measure towards embedded radio sources 
(e.g., Vogt & Ensshn 2005), spatially resolved ICM pres- 
sure fluctuations (Schuecker et al. 2004), and X-ray res- 
onant line scattering (e.g., Gastaldello & Molendi 2004). 
A direct observational link with our proposal that HBI- 
driven turbulence is present in cluster cores is via tur- 
bulent diffusion of metals throughout the core of the 
cluster. Rebusco et al. (2005) estimate the iron injec- 
tion rate into the Perseus cluster as a function of radius, 
and find that significant diffusion of the iron is required 
to explain the measured radial profile of the iron abun- 
dance. They estimate a diffusion constant within the 
core oiV ^ 2 y. 10^^ cm^ s~^ and show that this can be 
reasonably accommodated by a turbulent transport pic- 
ture. The source of turbulence is generally attributed 
to the central AGN, either by direct stirring of the ICM 
by AGN jets, or by AGN cosmic-ray driven convection 
(Chandran & Rasera 2007). Cluster-cluster mergers have 
also been considered as a source of turbulence (Fujita, 
Matsumoto & Wada 2004). 

A distinguishing feature of our model is that it makes 
a direct connection between the temperature gradient 
found in cluster cooling cores and turbulent diffusion. 
The inferred Rebusco et al. (2005) diffusion constant 
is readily accommodated by eq. (jlSp if ^ is order O.lr. 
Moreover, we also predict efficient turbulent transport of 
metals within the cooling core, but significantly reduced 
metal transport in the more isothermal regions of the 
cluster. Thus, the fact that the high metallicity cores 
appear to be limited to the region where a temperature 
gradient is maintained (e.g., Baldi et al. 2007; Snowden 
et al. 2008) is neatly accounted for. 

This scenario would require an associated convective 
heat flux larger than is seen in numerical simulations of 
noncooling HBI systems (Parrish & Quataert 2008). Es- 
timates of the ratio of the convective to local conductive 
fluxes based on equation (fT2|) and the above T) value are 
of order 10~^ to 10~^ for cluster cores. The noncooling 
HBI simulations, by contrast, yield values of order 10"'^. 
This will be an important point to reconcile in numerical 
simulations of cooling HBI flows. 

Thus far, we have focused exclusively on hot clusters 
where unbridled Spitzer conductivity is a more than suf- 
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ficicnt heat source for the cores. How does our picture 
relate to clusters more generally? Best et al. (2007) have 
studied a large sample of clusters from the Sloan Digital 
Sky Survey (SDSS) and gauged the role of radio-mode 
feedback in these systems under the very reasonable as- 
sumption that the relationship between radio-luminosity 
and mechanical power is not a strong function of envi- 
ronment. They find that, while embedded radio galaxies 
can more than offset radiative cooling in galaxy groups 
and small clusters, the AGN feedback hypothesis has sig- 
nificant energy problems when one considers large/rich 
clusters. They further show that (when parameterized in 
terms of the virial velocity of the cluster) the AGN feed- 
back hypothesis becomes untenable just about where the 
clusters are hot enough for thermal conduction to play a 
significant role. This coincidence is striking, and strong 
circumstantial evidence that thermal conduction plays 
an important role in regulating these hot clusters. 

5. SUMMARY 

Cooling flows in galaxy clusters show remarkably regu- 
lar temperature profiles. Naive cooling models, in which 
radiative losses are replenished by advected heat infiow 
(i.e., classical cooling flow models), fail badly, predicting 
an excessively large X-ray surface brightness and large 
amounts of cold gas (or powerfifl star formation). Di- 
rect heating from a central AGN is in principle a potent 
energy source, but why this should result in regular tem- 
perature profiles and temperature floors that scale with 
the virial temperature, is unclear. 

Temperature gradients in the hot gas would produce a 
large inward heat flux if Coulomb collisions alone were 
responsible for heat transport. This process, which has 
in fact been extensively studied, does not currently en- 
joy widespread support as an important generic heating 
mechanism. Heating by thermal conduction, when im- 
portant, seems to be overwhelming, and even if the match 
is good, thermal stability and fine-tuning remains an is- 
sue. For these reasons, researchers have often looked to 
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external sources of heat (e.g., the AGN) as the resolution 
of the cooling flow problem. 

We argue that the recent finding of Q08, that inward 
directing heat flow in a dilute plasma is unstable and 
strongly self-regulating, shows great promise as a mech- 
anism to maintain regular temperature proflles. In this 
work, we have proposed that the combination of a re- 
verse convective heat flux and field-line insulation main- 
tains critical stability conditions corresponding to nearly 
static stratified atmospheres. In essence, cluster cores are 
in a state of marginal stability, not to the HBI, but to the 
radiative HBI. Departures from this state should be sta- 
ble: an increased diffusive flux would produce a stronger 
reverse convective flux, and vice- versa. A fully insulated 
cooling core would draw field lines inward, opening a 
thermal path. As a by-product of convection, we expect 
efficient diffusion of metals in regions where a strong tem- 
perature gradient exists. This agrees with observations 
of enhanced ICM metallicities over precisely the same re- 
gion where the temperature gradient is maintained, but 
the associated convective heat flux would be larger than 
those inferred from nonradiative numerical simulations. 
Radiative HBI simulations are clearly needed. 

The basic scenario outlined in this paper — that cluster 
cores are in a state of marginal stability regulated by 
the radiative HBI can be tested numerically by direct 
simulation. In a subsequent paper, we will report the 
results of such a study. 
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